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' ABSTRACT 
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Aims. We present a new sample of active galactic nuclei (AGNs) identified using the 
catalog of the AKARI Mid-infrared (MIR) All-Sky Survey. Our MIR search has the 
advantage of detecting AGNs that are obscured at optical wavelengths by extinction. 
Methods. We first selected AKARI 9/im excess sources with F(9/J,m) /F(Ks) > 2 where 
Ks magnitudes were taken from the Two Micron All Sky Survey. We then obtained 
follow-up near- infrared spectroscopy with the AKARI /IRC to confirm that the excess 

5—i < 

is caused by hot dust. We also obtained optical spectroscopy with the Kast Double 
Spectrograph on the Shane 3-m telescope at Lick Observatory. 

Results. On the basis of these observations, we detected hot dust with a characteristic 
temperature of > 500K in two luminous infrared galaxies. The hot dust is suspected 
to be associated with AGNs that exhibit their nonstellar activity not in the optical, 
but in the near- and mid-infrared bands, i.e., they harbor buried AGNs. The host 
galaxy stellar masses of ~ 4 — 6 x 10 9 Mq are small compared with the hosts in 
optically-selected AGN populations. These objects were missed by previous surveys, 
demonstrating the power of the AKARI MIR All-Sky Survey to widen AGN searches 
to include more heavily obscured objects. The existence of multiple dusty star clusters 
with massive stars cannot be completely ruled out with our current data. 
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1. Introduction 



Many observations have found evidence of the presence of a large number of heavily 
obscured active galactic nuclei (AGNs). It is found that a si gnificant number of ultra- 
luminous infrare d galaxies contain AG Ns at their centers (e.g. 



S anders fe Mirabel 



1996 



Lutz et al 



1998 



Imanishi 



tions ( Maiolino et al 



et al 



1998 



2008). V arious hard X-ray and soft gamma-ray observa- 



Risaliti et al 



1999: 



Malizia et al 



2009 ) indicate that about 



80 percent of the AGNs in the local Universe are obscured. AGN synthesis models of the 
X-ray background post ulate their existenc e in order to explain the flat spectrum of the hard 
X-ray background (e.g. Ueda et al. . 2003|). T hese black holes contribute to the local black 



hole mass density (jFabian fc Iwasawa . 



1999), and are potentially important contributors 



to the growth of supermassive black holes throughout the history of the universe. However, 
the nature of this population, even in the local universe, is only poorly understood, because 
of the strong selection bias against finding them at optical wavelengths. 

Mid-infrared(MIR) AGN searches can overcome this obstacle by penetrating through 
dust extinction to identify most of the AGN population, includ ing Type 2 Seyferts 
and buried AGNs . The original IRAS 12/im active galaxy samples (|Spinoglio fc Malkanl . 



1989: 



Rush et al 



1993) provide an unbiased sample of local active galaxies. Using the 



ISOC AM parallel mode survey of 10 square degrees at 6.7^mfLW2 uLeipski et al 



(2005 



2007) succeeded in finding redder Type 1 AGNs as well as Type 2's. Several searches have 
also been performed u s ing the near- and mid-infrared bands in the Spitz er Space Telescope 
(e-g- 



Lacv et al 



2004; 



Alonso-Herrero et al 



2006 



Pollctta et al 



2006). 



AKARI performed an all-sky survey at 9 and 18/um as well as at four far-infrared (FIR) 
bands (65, 90, 140, and 160/zm). It provided improvements of about one order of magnitude 
compared to that of IRAS in both spatial resol ution and sensitivity in the mid-infrared 



Ishihara et al 



(2010). 



bands. The details of the survey are described in 

In this paper, we present the first results of our search for AGNs based on this AKARI 
MIR All-Sky Survey. We discovered two galaxies LEDA 84274 and IRAS 01250+2832, 
which have a compact hot >500 K dust component. The hot dust component may be 
heated by the central engine of the AGN, even though their optical spectra do not show 
any AGN characteristics. 

The observations, data reduction, and results are described in Section [2] In Section [3l 
we describe the multiwavelength properties of the two galaxies and in Section U the hot 
dust components we found in the galaxies are discussed. The discussion is presented in 
Section [SJ and a summary is given in Section [5J Throughout the paper, we assume a flat 
cosmology with fl — 0.3, A = 0.7, and H a = 70 km s _1 Mpc -1 . 



2. Observations and results 



The initial identification of the AKARI MIR All-Sky Survey sources involve s association 
with the Two Micron All Sky Survey (2MASS) catalog (jSkrutskie et all 120061 ). This search 
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highlights unusually red AKARI MIR sources with F (9 //m) / F (K s) > 2, at high Galactic 
latitudes, \b\ > 30° after excluding regions around the Large and Small Magellanic Clouds. 
To examine the origin of the excess of F(9/jim)/F(Ks) > 2, we performed follow-up obser- 
vations with the AKARI near-infrared(NIR) spectrometer. 

Using the AKARI spectra of the 2.5-5 fim wavelength range, AGNs can be distin- 
guished by their red continuum emission, while strong polycyclic aromatic hydrocarbons 
(PAH) emission is detected in star-forming galaxies. To measure the redshift and search 
optically for AGN or star-formation signatures, optical spectra were also taken with the 
Share 3m telescope at the Lick Observatory. 

During these follow-up observations, we discovered two buried AGN candidates that 
have steep red NIR continuum from hot dust, but do not show any AGN features in op- 
tical spectra, such as strong high-ionization emission lines. Table [T] summarizes the prop- 
erties of our targets, including the fluxes in the AKA RI All-Sky MIR and FIR Survey 



catalogs (jlshihara et al 



2010: 



Yamamura et al 



2009). LEDA 84274 is identified as the 



FIR s ource, IRAS 14416+6 618. Previous optical spectroscopy of LEDA 84274 (|Kim et al 



1995 



Veilleux et al 



1995) reported a low-ionization spectrum with strong Balmer lines, 



indicating an HII-region-like or star-forming galaxy at z=0.0377. Thus, the gas emitting 
these optical lines is photoionized by early-type stars. For IRAS 01250+2832, this is its 
first reported identification as a galaxy at z = 0.043. 



2.1. AKARI near-infrared spectroscopy 



We performed NIR spectroscopy using the InfraRe d Camera (IRC) 



Ohvama et al 



20071 ) onboard the AKARI satellite ( Murakami et al 



Ona 


<& et al. 




2007 


2007 


). We used the 



IRC channel NIR, which uses a 512 x 412 InSb array and the astronomical observation 
template Z4 (AOT Z4) designed for spectroscopy. AOT Z4 replaces the imaging filters by 
the transmission-type dispersers on the filter wheel to take NIR spectra. The NIR grism 
(NG) was set to cover the wavelengths of 2.5 - 5/im at a resolution of about 100. The 1' x 
1' slit was used to avoid any overlap with other sources. 

AKARI NIR spectra were obtained during December 2008 - July 2009. The total ex- 
posure times for both objects are 792 sec. We observed both objects twice in the AKARI 
warm (Phase3) mission period to correct hot pixels caused by the relatively high temper- 
ature (~ 40K) of the detector. 

The data were processed usi ng the IRC P h ase3 dedicated data reduction package, 
IRC_SPECRED Ver. 20090211 (johvama et all 120071) . Dark subtraction, linearity cor- 
rection, flatfield correction, and various image anomaly corrections were first performed. 
Multiple exposures were then coadded. After performing wavelength- and flux-calibration, 
the spectra of the object were extracted. The aperture size for the spectrum extraction was 
set to be 3 pixels (4.5"). Aperture corrections were also applied at the end of the process- 
ing. The present data reduction gives an uncertainty in wavelength of ~ 0.01/im and the 
error in the absolute flux is smaller than 20 percent. 

Figure U shows the calibrated spectra of LEDA 84274 and IRAS 01250+2832. Both 
objects display steep red continuum. We find emission lines of PAH at 3.3/im and Bra 
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Table 1. Data for LEDA 84274 and IRAS 01250+2832 



Property 


LEDA 84274 


IRAS 01250+2832 


R.A.(J2000) a 


14 42 34.88 


01 27 53.95 


Dec.(J2000) a 


+66 06 04.3 


+28 47 51.0 


Redshift 


0.0377 


0.043 


J(mag) 6 


13.756 ± 0.054 


14.774 ± 0.117 


if(mag) 6 


13.293 ± 0.090 


13.941 ± 0.131 


Ks(mag) b 


12.744 ± 0.083 


13.462 ± 0.169 


IRAS 12^m(Jy) c 


0.10 + 0.03 (2) 


< 0.11 (1) 


IRAS 25fim(Jy) c 


0.56 ± 0.03(3) 


0.28 + 0.05(2) 


IRAS 60^m(Jy) c 


2.19 + 0.20(3) 


0.52 ± 0.03(3) 


IRAS 100^m(Jy) c 


1.80 + 0.14(2) 


< 0.76 (1) 


AKARI 9/im(Jy) d 


0.084 ± 0.007 


0.105 + 0.010 


AKARI 18/im(Jy) d 


0.344 + 0.016 


0.182 + 0.017 


AKARI 65/im(Jy) e 


1.89 + 0.22(1) 


< 3.2* 


AKARI 90/im(Jy) e 


1.51 + 0.04(3) 


< 0.55 ,; 


AKARI 140/im(Jy) e 


1.79 + 0.26(1) 


< 3.8* 


AKARI 160/im(Jy) e 


< 7.5 4 


< 7.5* 


Lir(L )' 


1.8 x 10 11 


1.0 x 10 11 


Lir(L )* 


1.8 x 10 11 


0.9 x 10 11 


Ai(4000) h 


1.10 + 0.02 


1.63 + 0.08 



a The 2MASS coordinates. Units of right ascension are hours, minutes, and seconds, and units 

of declination are degrees, arcminutes, and arcseconds. 
b The 2MASS magnitudes in All-Sky Ext ended Source Catal og 

c Fluxes are from the IRAS Faint Source (jMoshir et al.l . ll992l l. Numbers in parentheses indicate 
the flux qualities. 



AKARI /IRC Mid- infrared All-Sky Survey Point Source Catalo g Ver. 1 (Ishihara 



et al 



2010) 



AKARI /FIS FIR All-Sky Survey Bright Source Catalog Ver. 1 jYamamura et all , 120091 ). The 
number in parentheses indicates the flux qualities. (3) denotes a clear detection, and (1) indi- 
cates the flux at the position. 

The 8-1000 /im infrared luminosity Lm, = AivD^Fm is computed usin g IRAS fluxes, F xa = 
(1.8 x 10" 14 )(13.56 x F12 + 5.26 x F 25 +2.54 x F 60 + 1.0 x Fioo) W m' ^Kim et allll995l ). We 
used IRAS upper limits when IRAS did not detect a source. 

The 8-1000 /im infrared luminosity is derived from the blackbody fits. See Sections 13. II and 13. 21 
for details. 

D n (4000) is the di scontinuity of the spect rum around 4000A. See Section \2. 2 1 for details. 



Detection limits in 



Yamamura et al 



(2010), which are estimated from the peaks of log iV — log S 



plots, corresponding to ~ 90 percent completeness. 



at 4.05/im for LEDA 84274. The measured line fluxes are given in Table [2] For IRAS 
01250+2832, the R-branch of CO ro-vibrational absorption is detected at 4.75/im. Upper 
limits to the fluxes of the PAH bands at 3.3/Ltni and Bra at 4.05/im are given in Table 2. 
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(a) LEDA 84274 
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(b) IRAS 01250 + 2832 
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Fig. 1. AKARI NIR spectra of (a) LEDA 84274 and (b) IRAS 01250+2832 in the observed 
wavelength frame. The thin vertical lines represent la errors. 

Table 2. AKARI NIR emission features of LEDA 84274 and IRAS 01250+2832 



Name 


Obs. Wave." Flux 6 


Obs. EW C 




(/j,m) (10~ 16 ergs s _1 cm -2 ) 


(Mm) 


LEDA 84274 


PAH 3.3fim 


3.416+0.010 437 ± 30 


0.018 ± 0.001 


Bra 4.05/^m 


4.200+0.010 178 ± 26 


0.060 ± 0.009 


IRAS 01250+2832 


PAH 3.3/im 


< 75 d 


< 0.003 d 


Bra 4.05pm 


< 100 d 


< 0.001 d 



a Observed wavelength. 

b Line flux and la error in the measurement. 
c Observed equivalent width in fj,m. 

Upper limits(3 a) are estimated with the line widths for LEDA 84274. 



2.2. Lick optical spectroscopy 

We performed optical spectroscopy of LEDA 84274 and IRAS 01250+2832 using the Kast 
Double Spectrograph on the Shane 3-m telescope at the Lick Observatory in 2009 March 
and August, respectively. This spectrograph has two separate parallel channels - one opti- 
mized for the blue and the other for the red. The 600/4310 grism and the 600/7500 grating 
were used in the blue and red spectrometers, respectively The wavelength coverage was 
3700 - 8300 A with a small gap at around 5500 A produced by the dichroic beamsplitter. 
The exposure times for both galaxies were 1500 sec. 

The spectra of both LEDA 84274 and IRAS 01250+2832 are shown in Figure [U We 
identify the strong emission lines in LEDA 84274 and IRAS 01250+2832 as shown in Table 
[3J On the basis of these identifications, we determine the redshift of z =0.0377 for LEDA 
84274, which is consistent with the redshift previously determined by iKim et al 



The line ratios of the present results are consistent with previous ones, while the emission 
line flux ratios relative to those of the Ha line in bluer part of the spectra are slightly 
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(a) LEDA 84274 
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Fig. 2. Optical spectra of (a) LEDA 84274 and (b) IRAS 01250+2832. In (b), the positions 
of the H/34861A and [O III] 5007A lines are indicated. 



larger than those in 



Kim et al 



(|1995() . The spectrum of IRAS 01250+2832 shows a series 
of absorption lines in the blue part of its spectrum. We identify these absorption lines as 
Call H+K, RS, Cal g-band, Hry, and Fe I as shown in Figure^b). Although Ha 6563A and 
[N II] 6548,6583A emission lines are detected, no other emission lines are detected. We use 
the absorption lines to estimate the redshift of z — 0.043. 

To examine the galaxy stellar population, we measure the 4000A discontinuity. This 
discontinuity is originally defined by 



L»(4000) 



(A 2 



^T) I\+ Fvd\ 



(A 2 + -A+) D Fvd\ 



(1) 



where (AT/, AJ, A+, A+) = (3750,3950,4050,4250) A ljBruzuaj Il983h . A definition using 
narrower con t inuum bands (AT, A^, A+, A^) = (3850,3950,4000,4100) was introduced by 



Baloghetal 



1999) and is widely used. The advantage of the narrower bands is that the 
index is less sensitive to reddening. We adopt the narrow-band definition, and denote this 
index as L> n (4000). The resultant D n (4000) indices are 1.10+0.02 and 1.63+0.08 for LEDA 
84274 and IRAS 01250+2832, respectively. On the basis of these indices, LEDA 84274 has 
the spectrum of a late-type galaxy, and IRAS 01250+2832 that of an elliptical galaxy. 
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7 



Name 


Obs. Wave." 


Flux 


Obs. EW b 


ratio with Ha c 


ratio with Ha d 




(A) 


(10 -16 ergs s" 1 cm -2 ) 


(A) 


This work 


Kim et al. (1995) 


LEDA 84274 


[U11J 


OO i u 


Qi 7 n+ 1 Q 




'iS 

U.JO 


"^O 








7 /i_i_n o 


08 




tip 




1 90 ^-1- 1 1 
IZU.Ozn 1.1 


on fi-L-f\ 9 


99 


1 8 
U. lo 


fOTTTl 




oi.yzn i.i 




1 K 


1 9 
U. 1Z 


fOTTTl 




94D 04- 1 1 

Z^U.UIL 1.1 






U.O 1 


Hel 


6094 


18.8+ 0.5 


3.2+0.1 


0.03 




[01] 


6536 


23.0+ 0.4 


4.4+0.1 


0.04 


0.04 


[01] 


6603 


8.5+ 0.4 


1.6+0.1 


0.02 




[Nil] 


6794 


61.1+ 0.6 


11.8+0.1 


0.11 


0.10 


Ha 


6808 


549.0+ 0.5 


106.1+0.1 


1.00 


1.00 


[Nil] 


6830 


178.0+ 0.5 


34.6+0.1 


0.32 


0.31 


[SII] 


6968 


78.9+ 0.4 


15.5+0.1 


0.14 


0.14 


[SII] 


6983 


65.1+ 0.4 


12.8+0.1 


0.12 


0.11 


IRAS 01250+2832 


[Nil] 


6828 


4.3+ 0.8 


1.5+0.3 


1.02 




Ha 


6844 


4.2+ 0.8 


1.4+0.3 


1.00 




[Nil] 


6864 


12.1+ 0.8 


4.2+0.3 


2.88 





Observed wavelength. A typical error of measurements is less than lA. 
Observed equivalent width. 

Emission line ratio with Ha in this work. A typ ical measurement error is less than 0.01. 



Emission line ratio with Ha in 



Kim et al 



( 19951 ). 



3. Multiwavelength properties 

The AKARI NIR spectroscopic observations display a steep red continuum in both galax- 
ies, while they have different characteristics in the optical and FIR. The following subsec- 
tions describe the properties of each galaxy. 



3.1. LEDA 84274 



The various properties of LEDA 84274 indicate that it currently has a high star- 
formation rate. The optical emission line ratios such as [NII]6583A/Ha vs. [OIII]5007A/H/3, 
[SII]6716A+673lA/Ha vs. [OIII]5007A/H/3, an d [OI]6300A/Ha vs. [O IH]5007A/H/3 are in 
the range of a star- forming gala xy according to 



Veilleux et al 



trum taken by 



Kim et al 



(|l995[ ). who used the spec- 



( 199,4) . detected at 3.416 is known to 

be associated with star-formation activity. FIR detections with IRAS and AKARI are also 
consistent with a high star- formatio n rate. The 20cm flux density in the NRAO VLA Sky 
Survey (NVSS: 



Condon et al 



19981 ) is reported to be 4.2±0.4mJy. The logarithmic ratio of 
the FIR to radio continuum flux densities log [F fir /(3.75 x 10 12 Hz)// l/ (20cm)] = 2.7 + 0.1, 
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where Ffi r is the flux between 42.5 and 122.5^nJ 1 l( Helou et al. . 1985 ). is consistent with 



2.4 ± 0.2 of the star-formation activity in n ormal, infrared, and lumin ous infrared galaxies 
including ultra-luminous infrared galaxies (panders fc Mirabel 119961 ). This relative weak- 
ness of radio emission in LED A 84274 is cl early inconsistent with the ratio (2.1) observed 
in local Seyfert galaxies (jRush et all 1996). 

Before we estimate the star-formation rate in LEDA 84274, we estimate the amount 
of reddening in the emission region, based on the Balmer line ratio, Ha/H/3. The equiva- 
lent widths of the Balmer lines are so large that we did not apply any correction for the 
underlying stellar absorption lines. The result is Ay = 1.3 mag , when we assume case 
B recombination with T e = 1 0000K (jOsterbrock fc Ferlandl 12006!) and the extinction law 
20011 ). while we get Ay = 3.5 mag using the recombination line 



( Wcingar tner fc Draine 



ratio of Bra/Ha. The difference can be explained by the deeper penetration of Bra pho- 
tons, suggesting that only the surface of the emitting region contributes to optical emission 
lines. 

We calculate the SFR using the PAH 3.3/xm band, which is not expected to be affected 
by an AGN. With the empirical ratio, £pah3.3/£/_r ~2-3x 10 -4 for starburst galaxie s 



(Mourietal 



1990T) and the conversion of the IR luminosity into SFR (jKennicuttl . 119981) 



the SFR of LEDA 84274 is estimated to be 20 - 30 M Q year" 1 . We estimate SFRs of 7 and 
25 M Q year" 1 using the optical Ha emission line with Ay = 1.3 and the NIR Bra emission 
line with Ay = 3.5, respectively. SFR(Bra) is consistent with SFR(PAH), while SFR(Ha) 
is lower than the others, which can be explained similarly by a highly non- uniform Ay. 
The infrared luminosity calculated from the IRAS measurement is Lir = 1.8 x 10 11 Lq, 
which corresponds to SFR(IR)=31 M year -1 . 

The spectral energy distribution (SED) of LEDA 84274 is presented in Figure El A 
SEP model is also o verplotted in Figure |3][a) . For LEDA 84274, the Sc galaxy template 



(Polletta et al 



20071) agrees with the photometric points in the J— , H-, and if 5— bands. 



To fit the AKARI NIR spectrum, we assume that the spectrum is a summation of the 
Sc galaxy template and dust with a single-temperature modified blackbody. The model is 
defined as 



T\ C\ Tgalaxy + C 2 B\(T dust )(l - exp(-rx)), 



(2) 



where -Fgaiaxy is the Sc galaxy template, -Ba (Tdust) is the Planck function with dust tem- 
perature Tdust and t\ is the dust optical depth. For the frequency dependence of the dust 
optical depth, we adopt t\ = CgA - ^. The parameters C\, C2, Tdust, C3, and j3 are derived 
by x 2 fitting the AKARI NIR spectrum. The fit indicates that 530K < T dust < 670K if 
optically thick, or 470K < T dust < 580K if optically thin with n Mm < 1 and ,0=1—2. 
Figure [3] (b) shows a close-up of a fit with an optically thick blackbody of Td us t = 600TT. 

In addition, we must add a cool component of dust with Tdust = 93 ± 5 K and — 
to explain the IRAS and AKARI photometric data at wavelengths longer than 12/xm. For 
AKARI 9/im and IRAS 12/im data, we need to add one more component of temperature 
Tdust ~ 300K. There is no strong constraint on either the temperature or the emissivity. 



1 Ff ir = 1.26 x 1CT 14 x [2.58/„(60Aim) + /„(100/im)] where /„ are the flux densities in Jy, and 
Ffi r is in Wm" 2 . 
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The fitting results are shown in Figure[3](a). From the blackbody fits, we derived a total 
infrared luminosity of L\r = 1.8 x lO 1 1 !^ between 8 an d 1000 /im , which is consistent with 
the Lir derived using the method of 
is dominated by the cool component. 



Kirn ct al 



1995). We note that the total luminosity 



3.2. IRAS 01250+2832 



The SED of IRAS 01250+2832 is shown in Figure H It is difficult to explain the fiat 
SED a t wavelengths bet ween 9/im and 60/im with normal spiral and elliptical galaxy 



models( Silva et al 



19981 ). The spiral and elliptical galaxy templates do not have strong 



enough infrared emission. It is difficult to explain the global SED with starburst galaxy 
templates such as M82 and Arp 220 that have a strong excess in the FI R. Here, we assume 



Silva et al 



1998)) 



that the host galaxy is dominated by an old stellar population (13Gyr; 
based on the strong £) n (4000) index. 

When we also fit the AKARI NIR spectrum using Equation [5] and the elliptical galaxy 
template, we need to introduce a blackbody component of temperature of 510 ± 20K in 
optically thick conditions as shown in Figure[H(b). Another blackbody component of 100K 
is needed to explain the IRAS 60^m flux. The higher fluxes of the model relative to those 
observed at 9 and 12 /jm can be attributed to 9.7/im silicate absorption, which is not 
included in this model. However, it is expected to be present based on the CO absorption 
measured at 4.75/mi. The conclusion that the SED needs two blackbody components does 
not change much even if we use spiral galaxy templates. The infrared luminosity Lir is 
derived to be 9.1 x 10 10 Lq, which is consistent with 1.0 x 10 11 Lq using the method of 



Kim et al 



(1995). We note that the contribution at wavelengths longer than 60/im to the 
total luminosity is expected to be small for IRAS 01250+2832. 

While there are no signs of H/3 and [OIII] emission lines in Figure HJb), Ha and [Nil] 
lines are clearly detected. The [NII]6584A emission line is stronger than Ha, suggesting 
that IRAS 01250+2832 is a low-ionization nuclear emission-line region (LINER) galaxy. 
However, the other indicators of LINERs, [OI]6300A and [SII]6716A+673lA emission lines, 
are not seen because of the weakness of the lines, relative to the stellar continuum. The 
3cr upper limit to the H/3 flux is estimated to be 3.0xl0 -16 ergs s _1 cm~ 2 . This upper 
limit is not useful to estimate the reddening from the Balmer ratio. Without applying 
any reddening correction, we calculate an Ha line luminosity of 4.7xl0 5 Lq, which is 
comparable to that of low -luminosity AGNs (LLAGNs) fou nd in the optical spectroscopic 
survey of nearby galaxies (|Ho. Filippenko. fc Sargentl . 119971 ). A comparison with LLAGNs 
is discussed in Section [5] 



4. Hot dust with T dust ~ 500 - 600K 

As shown in Section [3J both galaxies show near-infrared continuum from hot dust with 
Tdust = 500 — 600K, although they differ in other ways. The origin of the energy source 
heating this dust, however, remains unclear. 

Fits indicate that the luminosity of this hot dust component is 2.7 x 10 10 L Q and 
6.6 x 10 10 L for LEDA 84274 and IRAS 01250+2832, respectively. If we assume one 
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(a) 



10 



10' 



00 




Fig. 3. Spectral energy distribution of (a) LEDA 84274. The blue line shows the AKARI 
NIR spectroscopy. Squares and triangles represent the AKARI All-Sky Survey and IRAS, 
respectively. Crosses s how the 2MASS phot ometry. The solid line(black) is a model tem- 



plate of a Sc galaxy (jPolletta et al 



2007 ). The model line normalized at the 2MASS 



photometry. A dashed line (red) represents the Sc galaxy template plus three blackbody 
components, the temperatures of which are 600K, 300K, and 93K. Dotted lines represent 
each blackbody component, (b) Close-up of the 2-5 /im range of LEDA 84274, to clearly 
show the model and the blackbody component with 600K. See text for details. 

(b) 





Rest wovelength( J u,m) 



Rest wavelength(/arn) 



Fig. 4. Same as Figure E] except for IRAS 01250+2832. A model template of an elliptical 
galaxy (jSilva et al.1 . 119981 ) is used. Two blackbodies in this figure have temperatures of 
510K and 100K. 
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single optically thick spherical emitting region, the size of the emitting regions would be 
0.7 pc and 1.6 pc in diameter for LEDA 84274 and IRAS 01250+2832, respectively. One 
possibility is that we are observing an OB star cluster that has a huge luminosity of 
-Lir > 10 Lq in a compact volume of a cubic parsec, and that heats up the surrounding 
dust to >500K. However, the relaxation time of the dynamically bound system can be 
estimated to be 10 6 — 10 7 years. Thus, this situation is unlikely to be observed. 

A more likely possibility is that the object is a compact central supermassive black hole 
surrounded by a hot accretion disk. However, no optical spectrum shows any evidence of 
the hard photons produced by an AGN. The various emission -line ratios of LE DA 84274 in 

( 19951 ) concluded. 



Kim et al 



the optical band show that it is an HH-region-like galaxy as 
IRAS 01250+2832 shows no emission lines except for Ha and [Nil]. The weak- and non- 
detection of PAH 3.3/im emission may support the presence of an obscured AGN. The 
equivalent widths of the PAH 3.3 /mi emission feature are 0.018 /jm and < 0.003 fim for 
LEDA 84274 and IRAS 01250+2832, respectively. Observations of starbursts show equiva- 
lent wid ths of PAH 3.3 lira with an average valu e of 0.1/zm, and they are never smaller than 



1986 



Imanishi ct al 



20081 ). Thus the PAH equivalent width is diluted 



0.04um (|Moorwood , 

by strong featureless NIR continuum in IRAS 01250+2832 and completely overwhelmed in 
LEDA 84274. 

The excess featureless near-IR continuum may alternatively represent multiple dusty 
OB star clusters, in whic h extraordinar i ly intense starbursts are occur r ing. This is s een in 



some interacting galaxies (| Joseph et al 



19841: Uoseph fc Wright! Il985l) . iHunt et alJ (120021) 

also argue that the hot dust that produces the NIR color excess must be heated by massive 
stars. In this case, the weak PAH emission can be produced by the ha rder radiation fields 
of massive star s in these environments that can destroy PAHs (e.g. iBeirao et al 



2006 



Wu et al 



20061) . From Spitzer imaging data of nearby irregular galaxies, excess emission 



at 4.5 /im from hot dust was found mostly in high-surfac e-brightness HII reg ions, imply- 
ing that massive stars are the primary source of heating (jHunter et all 120061 ) . For IRAS 
01250+2832, multiple dusty star clusters are less likely to be present, because most of dusty 
star clusters are surrounded by diffuse emission of PAHs as well as optical emission lines. 
IRAS 01250+2832 does not show any evidence for such activities. 

If the assumption that the dust emission is optically thick is not correct, the NIR red 
continua of LEDA 84274 can be explained by a dust temperature of ~ 480 — 580K and 
emissivity of with j3 = 1 — 2. In spiral galaxies that normally display the thermal 

(|2003[ ) detected non-stellar NIR excess continuum 



Lu et al 



emission of optically thin dust 
with a temperature of ~ 1000K(/3 = 2) and concluded that the NIR excess continuum 
originates in the interstellar matter of the galaxies based on the linear correlation between 
emission from aromatic carbon and the excess. In this case, the NIR excess has a luminosity 
of only a few percent of the FIR luminosity. However, the ratio of the NIR luminosity to 
the FIR luminosity is 0.15 and 0.75 for LEDA 84274 and IRAS 01250+2832, respectively. 
Thus, the emission from the hot dust in LEDA 84274 and IRAS 01250+2832 seems to be 
different from those of spiral galaxies. 

We suspect that the hot dust associated with obscured AGNs. But the possibility 
of multiple dusty star clusters in both objects, especially LEDA 84274, cannot be ruled 
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out. Observations at other wavelengths can resolve the question. In particular, X-ray ob- 
servations can provide important information on the energy source, but both galaxies 
are not detected in the ROSAT all-sky surv ey faint source cat alogue (an upper limit of 
1.3 x 10 -13 ergs cm -2 in the 0.1-2.4 keV band; 



Voges et al 



2000) . Hard X-ray observations 



will thus be important. 



5. Discussion 



As discussed in Section the compact hot dust component that we found in LEDA 84274 
and IRAS 01250+2832 is likely to surround an AGN that is heating it, even though other 
possibilities cannot be completely ruled out. In this section, we discuss these dusty AGNs. 

The optical spectra display no strong emission lines (especially high-excitation ones) 
from narrow-line regions. Thus, no narrow-line region on a 10-1000pc scale is evident 
in these galaxies. IRAS 01250+2832 has a low-ionization spectrum that is typical of 
LLAGNs, as discussed in Section l3~2l The low-ionization state spectrum in LLAGNs may 
be produced by bla ck-hole ac cretion at a very low rate (with a bolometric luminosity of 
Lboio < 4 x 10 10 L (|hoI . 120081) 1. However, LEDA 84274 and IRAS 01250+2832 have much 
higher luminosities of Lir > 1 X 10 11 Lq at infrared wavelengths than LLAGNs. Large 
optical extinction of the central engine is a plausible explanation of their optical spectra. 
These findings indicate that the AGNs are surrounded by a large amount of dust, so that 
ionizing UV radiation from the AGN is blocked along virtually all lines of sight. In ultra- 
luminous and lumino us infrared galaxy cores, similar situations are sometimes seen (e.g. 



Imanishi et al 



20101). 



Since the observed radiation from AGNs is often reprocessed and re-emitted at infrared 
wavelengths, the intrinsic SED is altered substantially. Wide spectral coverage is necessary 
to de termine the Eddington ratio (-Lboio/^Edd, or equivalently iboio/AfBH) (e.g. [Peterson, 



1997T ). Although it is difficult to differentiate the AGN component from the overall SED, we 



estimate the black-hole mass by making the a simple assumption that the AGN luminosity 
is nearly equal to the luminosity from the 500-600K dust. Assuming spherical accretion 
toward the black hole with a radiative efficiency of 0.1, the supermassive black hole mass 
from the luminosities of the hot dust is at least 8 x 10 6 M and 2 x 10 7 M Q for LEDA 
84274 and IRAS 01250+2832, respectively. These do not differ significantly from the black 
hole masses estimated from the 9/xm luminosit i es usi ng the bolometric correction evaluated 
using a Type2 QSO template in 



Polletta et al 



(|2007[ ) 3 Thus, we adopt the black hole mass 
estimates from the luminosities of the hot dust in the following discussion. If they have a 
higher radiative efficiency, the masses will be lower. 

The masses of these host galaxies can be calculated to be 6 x 10 9 M Q and 4 x 10 9 M Q for 
LEDA 84274 and IRAS 01250+2832, respectiv ely, using the 2M ASS K s photometry and 

. The mass ratio of the 



the mass-to-luminosity ratio, Mhost/-^ic s ~ 1 (|Cole et al 



2001 



s 



The hot dust luminosities estimated from the 500-600K blackbody are 30 percent higher and 
4 percent lower than those estimated from the bolometric correction, for LEDA 84274 and IRAS 
1 250-L2232^j£saeciisiely. 



In 



Cole et al 



12001), mean stellar mass-to-light ratios are estimated to be 0.73 and 1.32 



Mq/Lq for the Kennicutt and the Salpeter initial mass functions, respectively. 
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black hole to the stellar component in the host galaxy is 1.3 x 1CP 3 and 5 x 10~ 3 for LEDA 
84274 and IRAS 01250+2832. The results are consistent with the local relation between 
the mass of the central black hole and the stellar mass of the surround ing spheroid or the 



2004). The mass of both 



bulge in nearby galaxies, M^jMbuige = 1.4 x 10 3 (|Haring fc Rix 
host galaxies is relatively low, i.e., 6 x 10 9 M Q and 4 x 10 9 M Q , for LEDA 84274 and IRAS 
01250+2832, respectively. Large samples of AGNs from the Sloan Digita l Sky Survey have 
host galaxy masses of logM ffo ; = 9.5 — 12 with D n (4000) = 1.2 — 2.2 (jKauffman et al 



2003). Our sample with Af host = 4 - 6 x 10 9 M and D„(4000) = 1.1 - 1.6 seems to be the 



least massive population that may harbor an AGN. 
6. Summary 

By combining data from the AKARI MIR All-Sky Survey and NIR spectroscopy, we have 
discovered two luminous infrared galaxies at 0.04 that have hot dust with a temperature 
of >500K. The hot dust is likely to be associated with AGNs, but the possibility that the 
hot dust is heated by multiple dusty star clusters of massive stars cannot be ruled out by 
our current data. If they are AGNs, their emission must be buried, because present AGN 
signatures are seen only in the NIR and MIR, not in the optical. If they have AGNs, the 
estimates of the black hole masses of ~ 8 — 20 x 10 7 M Q , and the host galaxy masses 
of ~ 4 — 6 x 1 9 Mq also indicate th at the population is less massive than optically 



selected AGNs ( Kauffman et al 



2003). These objects were missed in previous surveys, 
demonstrating the power of AKARI MIR All-Sky Surveys to obtain a more complete view 
of the entire AGN population. 

We have initiated a long-standing program to study AGNs detected from the AKARI 
MIR All-Sky Survey. Multiwavelength observations will provide larger samples of buried 
AGN candidates that are similar to those we present in this paper, examples of which we 
plan to present in the near future. 
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